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Abstract 

We investigate the integrability of Nonlinear Partial Differential Equations (NPDEs) . The concepts are devel- 
oped by firstly discussing the integrability of the KdV equation. We proceed by generalizing the ideas introduced 
for the KdV equation to other NPDEs. The method is based upon a linearization principle which can be applied 
on nonlinearities which have a polynomial form. We illustrate the potential of the method by finding solutions of 
the (coupled) nonlinear Schrodinger equation and the Manakov equation which play an important role in optical 
fiber communication. Finally, it is shown that the method can also be generalized to higher-dimensions. 
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1 Introduction 



The conditions under which Nonlinear Partial Differential Equations (NPDEs) can be solved are even in one dimension 
not well understood JjJ. Roughly speaking the majority of the integrable systems can be classified in three main groups. 
In the first of these groups are those equations which can be reduced to a quadrature through the existence of an 
adequate number of integrals of motion. In the second class are those equations which can be mapped into a linear 
system by applying a number of transformations (hereafter to be called C-integrable) Q . The last group consists of 
differential equations which can be solved by Inverse Scattering Transformations (1ST). In the following, we will call 
equations which can be solved by inverse scattering methods "S-integrable" . The discovery of the 1ST has lead to a 
considerable progress in understanding the topic of integrability since this technique made it possible to investigate 
the integrability of large classes of NPDEs systematically [|| . Recently, it has been discovered that NPDEs which can 
be solved by inverse scattering techniques can also be transformed into linear differential equations j| . We therefore 
have the feeling that NPDEs which can be solved by inverse scattering methods form a sub-class of the C-integrable 
NPDEs. 

Another important consideration is that almost all the work on the integrability of NPDEs has been carried out 
in one space dimension only. Although the inverse problem of the Schrodinger equation can be generalized to three 
dimensions, the method is far too complicated to solve higher dimensional NPDEs. An alternative is the 9-approach 
which is also successfully generalized to N-dimensions (see for instance Ref.||). Nevertheless, for both these methods 
the existence of the obtained solutions is difficult to prove. The concept of C-integrability however, has the potential 
to be generalized easily to dimensions higher than one. In this paper, we will demonstrate a simple method based upon 
linearization principles that generalizes the concepts of C-integrability and S-integrability and enables us to compute 
general solutions of large classes of NPDEs by solving a linear algebraic recursion relationship. The result suggests that 
this method can also be generalized to systematically derive large classes of solutions of higher-dimensional NPDEs. 

In this paper we aim to find large classes of integrable equations which can be solved by linearization. Since it is 
not clear what integrability exactly means, we use in this paper the heuristic definition that an NPDE is integrable if 
given a sufficiently general initial condition, we can find analytic expressions the time-evolution of the solution. For 
NPDEs which can be solved by inverse scattering techniques, this notion is equivalent with the existence of N-soliton 
solutions, since it is implicitly assumed that the obtained solution can be expanded on a Fourier basis Qj. In this 
paper, we will show that the condition of expansion in a Fourier can be replaced by an arbitrary other complete set 
of basis-functions. 

We present the following novel results. Firstly, we derive a simple method to find general solutions of large 
classes of NPDEs. Secondly, we show that the integrability of these NPDEs is guaranteed if the nonlinearity can 
be expanded in the same basis- functions as the linear part, and secondly, if the dispersion relationship associated 
with the linearized problem can be solved. The method is firstly illustrated by discussing the integrability of the 
KdV equation in Section 2. In Section 3, the concepts derived for the KdV equation are generalized to discuss the 
integrability of general NPDEs. Finally, in Section 4, the results are applied to investigate the integrability of the 
coupled nonlinear Schrodinger equation and the Manakov problem which play an important role in the field of optical 
fiber communications. The paper is concluded with a discussion. 



2 The integrability of the KdV equation 

In order to illustrate the machinery developed throughout this paper we firstly discus the integrability of the KdV 
equation as an example. The integrability of the KdV equation is a well-studied problem j|. This makes the KdV 
equation an ideal object to test the validity of newly developed ideas with respect to the integrability of NPDEs. We 
will introduce our methods on the integrability of NPDEs by discussing the existence of N-soliton solutions for the 
KdV equation which is given by: 

u t +u xxx = 6u x u (1) 
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We try to find solutions of Eq.([l|) by substitution of the following Fourier series: 

oo 

u(x, t) = J2 A n (? n{kx ~ ut) (2) 

n=l 

If we substitute the solution u(x,t) into Eq.(||), we obtain: 

oo oo n— 1 

(nuj + k 3 n 3 ) A n e m{kx - Ut) = -6k ^ ^ lAtA^e™^'^ (3) 

n=l n=l 1=1 

We can now determine the coefficients A n by deriving a recursion relationship. This can be achieved by comparing 
the exponential functions in Eq. ([}]). If we compare all the terms for which n = 1, we find: 

(uj + k 3 ) A^ 1 ^^ =0. (4) 

For a nonzero A\ , we find that Eq. (Q) is satisfied if: 

w = -k 3 (5) 
If we put n = 2, we can determine A2 by solving the following relationship: 

(2cj + 8k 3 ) A 2 e 2l{kx - Ut) = -GkAiAie 21 ^-^ (6) 
If we use the dispersion relationship (||), we find that A2 is given by: 

A 2 

m = -£ m 

By repeating this procedure, we can compute all the expansion coefficients A n of the solutions u(x,t). In general, all 
the coefficients A n can be computed by solving the following linear algebraic problem: 

L^(k)A n = R^{k) (8) 

The operators L^ n \k) and R n (k) in Eq.(||) are given by: 

n-l 

L {n) (k) =n[n 2 -l]k 3 R {n \k) = -Qk^lAiAr^i (9) 

1=1 

If we compute all the coefficients A n by using Eq.(||), we than obtain the Fourier expansion of u(x, t) for which the 
first terms are given by: 

A 2 '*A 3 

u(a:, t) = Are^-^ - * e^**-"*) + £fl e -3i(*x-<rf) + . . . / 10 ) 
k 4fc 4 



If we substitute k = 2i/3 and A\ = Ad/3 into Eq. ( [To| ) we find: 



P 



U(£B, i) = J (12) 



By carrying out the summation in Eq.(]ll|), we can formulate this equation more compactly: 

(l + | e -2(/3*-4/3^ 
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Hence, if we put: 



2 V ' V~c \ P 

we can simplify Eq.(|12|) one step further to: 



i(x, t) — — -sech 2 |--y/c(a; — ct + xo) | 



(14) 



Eq.([lJ) describes the well-known KdV soliton. 

What did we learn from this simple exercise? At first, the KdV equation has solutions because of the special 
structure of the nonlinearity. If we substitute the special solution (j^) in the nonlinear part of the KdV equation, we 
find that we can expand the nonlinearity in the same basis functions as the linear part: 



oo n— 1 

6u x u = O n e in ( fa - w " ; D n = -6k £ IMA^-i (15) 

n=l 1=1 

This guarantees that we can find an iteration relationship for the expansion coefficients A n . As we will see later, 
we do not have to restrict to a Fourier expansion of the solution only. In principle this method works for any set 
of basis- functions as long as we can expand the nonlinearity in the same basis- functions as the linear part. In the 
following, we will show that the structure of the nonlinearity of the KdV equation enables us to construct the Fourier 
expansion of the N-soliton of the KdV equation. In order to systematically solve these solutions it is illustrative to 
discus also the two-soliton solutions which are assumed to have the following series expansion: 

Z\—X 



u(x,t)= ]T C( m ,M 2 )e^ lfclzl+ ^ fe2 ^ 1 J { i 2) (16) 
-, \ z 2 — x j-^ 

If we substitute Eq.([l6|) into the KdV equation ([!]), we obtain the following result: 



fll.fl 2 — l 

oo fJ-i— I5M2 — 1 

I"l,/i2 = l »?1,»)2 = 1 



(17) 



where 



L(" 1 '" 2 )(fc 1 ,fc 2 ) = ^]n l [n 2 -l]fc l 3 M ( - n ^\k ll k 2 ) = Y. n ^ ( 18 ) 

i=l i=l 

We solve Eq.(|l7|) by comparing equal exponential powers on both sides. This can be done by defining a parameter 
r = fii + fi2 and subsequently comparing the powers for T = 1, 2, 3, • • • . We firstly discus the case in which T = 1 in 
which only the coefficients C(l, 0) and C(0, 1) contribute: 

[ui + k'f] C(l, 0)e lklZl + [«a + fc£] C(0, l)e lk2Z2 = (19) 

If we put C(l, 0) = A\ and C(0, 1) = A 2 , we find that the following linear dispersion relationships must be valid: 

w(fci) = -hf and uj(k 2 ) = -k\ (20) 



4 



Once the linear dispersion relationships are determined and if the coefficients C(1,0) and C(0, 1) have taken their 
values A\ and A 2 , we can compute all the other coefficients C(/z, rf) by applying the following linear recursion relation: 

L^^\k 1 ,k 2 )C(fi 1 ,fi 2 )^R^ 2 \k 1 ,k 2 ) (21) 

where 

R^^(k u k 2 ) = -6 Mim ' m \ k ^ k 2)C(pi~Vi^2-V2)C( m ,r l2 ) (22) 

Equation (21) has a similar structure as Eq.(|8|). In principle Eq.(^l|) provides an efficient tool to compute all the 
coefficients C(ii,t]). We can easily generalize this result to the N-soliton case by assuming that the solution u(x,t) 
takes the following form: 

r z . = x _ 

oo \ k\ 

u{x,t)= J2 C(^ 1 ---fi N )e z ^ lklZl+ - + ^ NkNZN ' ) I : (23) 

We can determine the non-zero coefficients C(/^i • • • /ijv) by substituting Eq.(p3[) into the KdV equation ([!]): 



(24) 



J2 L^ 1 " , "* r) (fci • • • k N )C(p! ■ ■ ■ flN y(p^^+-+^ k « z «) 

Mi -"Mjv — 1 

OO /J,i — 1- ■ -/^jv — 1 

= -6 ^ ^ M^--^)(fc 1 ...fc w )C( Atl -r, 1 ...^-r yjv )C(?7i---w)e^ 1 ' ClZl+ --- + ' 1Nfc,v2N) 

where 

JV AT 

... kN )=J2 «i[n? - l]fcf ; M(" 1 -"~)(fc 1 • • • fcjy) = ^ (25) 

i=l i=l 

If we use that = -kf, (i e 1 • • • N) and A x = C(l, 0, 0, • ■ • , 0), A 2 = (7(0, 1, 0, • • • , 0), • • • , A N = C(0, • • • , 0, 1), 

we find that the expansion coefficients of the N-soliton solution for the KdV equation can be computed by solving the 
following linear relationship: 

L(Mi-w)( fcl . . . k N )C(m ■■■(i N ) = R^-^ih ■••k rf ) (26) 

where 

Hi — l---/ujv — l 

R^-^\k!---k N ) = -6 M^-^\k 1 ---k N )C( M - m ---ii N -r lN )C{T ]1 --- m ) (27) 

■qx— W=l 

From the exercise performed in this section we can conclude that general solutions of the KdV equation can be obtained 
by solving Eqs.(|2^). This implies that the KdV equation can be transformed into a simple linear algebraic equation 
in the coefficient space. We can conclude that the KdV equation has N-soliton solutions because the following two 
conditions are satisfied: 

• The structure of the nonlincarity of the KdV equation guarantees that the equation has solutions of the form 
(PI). This result implies that the coefficients R^"^{k\ ■■■k N ) exist. 
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• L^-^iki ■■■k N ) is not equal to zero if ki ■ ■ ■ k N ^ and n x ■ ■ ■ n N ^ 0. This implies that L^ 1 "^ (ki ■ ■ ■ k N ) 
has an inverse. 

In the following section we will show that a similar condition must hold for more general nonlinear NPDEs. In the 
following section it is shown that the concepts derived for the KdV equation can be generalized to large classes of 
NPDEs. The results obtained in this section are derived by assuming that the solution of the KdV equation can be 
expanded in Fourier basis functions. In the following section, it will be shown that similar principles apply for general 
basis functions. 



3 Generalizations 



In this section we will present more general results with respect to the integrability of nonlinear evolution equations. 
This will be done by generalizing the results obtained for the KdV equation. In this section, we focus on NPDEs of 
the following type: 



C[u(x,t)] = Q[u(x,t)} 



(28) 



In Eq(^8|), the function u(x, t) is a M-component vector function having entries Ui(x, t). The operator C[ ■ ] is assumed 
to take the following form: 



C[u(x,t)} = 



il 



d_ 

dt 



n=l 



(n) 



gn 



u(x, t) 



(29) 



The matrices A^™) in Eq.(|29|) are MxM matrices and I is the identity matrix. As concluded from the previous section, 
integrability puts strong constraints on the nonlinearity represented by the operator Q. As a necessary condition for 
the integrability we require that if a solution of Eq. (J2q) has the following form: 



u(x,t) = ^ C(fii • ■ • fi N ) exp 



N M 



% ^^^^^ ^ ^-L'pJ^'f^Z'p 



r—l s = l 

than, the operator Q must satisfy the following property: 

oo 

Q[u(x,t)}= ■ • • Mw) eoqp 

fJ-l ■■■MiV = l 



U)(k rs ) 



N M 



\j ^ ^^^^^ f-^T" ^^7* S 



r—l s—1 



(30) 



(31) 



where C(/ii • • ■ (J-n) and R(/ii • • • /j-n) are M-dimensional vector functions. Similarly as for the KdV equation, the 
vector function R(/ii ■ ■ ■ (J-n) is specified by the nonlinearity. In other words: we require that given a solution of the 
form (|3C|), the nonlinear operator Q[u(x,t)] can be expanded in the same set of basis functions as C[u(x,t)]. In the 
previous section, we have shown that the nonlinearity of the KdV equation satisfies this condition. In general large 
classes of nonlinear operators will have the property ( |3l| ) and among them we are especially interested in the sub-class 
P which plays an important role in nonlinear optics: 



, , du du 



d p u 



dx^dtP-i 



(32) 



where Pjv are polynomials of order N. If we let act the the linear operator C onto the solution (|30|) we obtain the 
following relationship: 



N M 



K 



£[U(X, *)] = I E E VMkrs) + J2 A 



(n) 



r—l s—1 



n=l 



N M 



u(x,t) 



(33) 



G 



From this result, we can identify a matrix L^ 1 ' /iJV ' (fcy) which is given by: 



N M 



K 



N M 



(34) 



This result implies that the coefficients C{ni • • • /ijv) which determine the solution ( |30| ) can be determined by solving: 

L^-^)(%)C(/i 1 --- AtiV )=R(Mi---/iiv) (35) 

The coefficients C(l, 0, 0, • • • ,0), C(0, 1, 0, • ■ • , 0), • • • , C(0, • • • , 0, 1) are determined by the initial condition. 

In principle we expand the solution u(x, t) in an arbitrary set of basis-functions. Suppose as an example a function 
u(x,t) which can be expanded in the set of basis-functions f( n >(x,t\k,u>): 



(36) 



We define the set S as the basis-functions: 

S = t\k, ^J^ix, t\k, u), ( Xj t \k, u), ■ ■ ■ }, 

which have the following properties: 

I: if/W(x,t|fc s «) G S => jLfM(x,t\k,u) = & n (k,u)fM(x,t\k,uj) 

II: if/W(x,t|fc s «) G S => ^f( n \x,t\k,u) = p n {k,u)f {m Xx,t\k,uj) 
III : if/<»(x,i) G S and / (m) (x,i) G S => f in) (x,t)-f ( - m \x,t) G S. 

The properties I and II guarantee that C[u(x, t)] can be expanded in basis functions p n >(x, t\k,w): 



(f^ m \x,t) G S) 
(fW(x,t) G S) 



C[u(x,t)} =^£< n >a n /< n >(x s t|fc,w) 



(37) 



(38) 



(39) 



Where the precise structure of the operator ZA n ) is determined by the linear differential operator C. Property III 
in Eq.(|38|) guarantees nonlinearities of the type P can be expanded in the same basis functions f^ n \x,t\k,u>). If the 
nonlincarity represented by the operator P can also be expanded in the same basis-functions f^(x, t\k, u): 



Q[u(x,t)]=J2R n fW(x,t) 



n=l 



than, we can compute the expansion coefficients a n by solving the relationship: 



(40) 



(41) 



where a\ is determined by the initial condition. Of course, we can generalize this result further by replacing Eq. (J30|) 
by: 



N M 



(42) 
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The structure of the solutions proposed in Eq.([42|) is in fact a generalization of Eq.(^o|). If we replace (x, t\kij , u>ij) 
by exp[ifJtikijZij], the form ( |30) i s retained. Following a similar approach as in the case of Fourier basis functions, 

38l) hold for the solution (E2I), the linear part of the differential equation acts on the 



we find that if the conditions 
solution (fH) like: 



N M 



K 



1=1 j=i 



N M 
*=1 J=l 



u(x, t) 



(43) 



where it is assumed that dtf^(x,t\kij,0Jij) — &ijf( % '(x,t\kij,u)ij) and d x f^(x,t) = kijf^'{x,t\kij,u>ij). This rela- 
tionship enables us to identify an operator LW)(£>y, kij) according to: 



L<«>(w W ,Ay) 



z— 1 j=l n— 1 



N M 

i=l i=l 



(44) 



If we moreover assume that the operator Q is of the class P so that: 

00 N M 

Q[u(x,t)]= J2 M^f-^HHf^^Xhj^ij) 

than the expansion coefficients are determined by the following linear iteration series: 

L (y) (a>y-, kij)C(jj,i ■ ■ -fi N ) = K(pi ■ ■ ■ n N ) 



(45) 



(46) 

From this result we can conclude that we can transform Eq.(|28|) into Eq.(|46|). We can conclude that an NPDE of the 
form Eq. ( |28| ) is integrable if the following two conditions are satisfied: 

• The nonlinearity must have such a structure that it can be expanded in the same basis functions as the linear 
part. In other words, the nonlinearity must guarantee that Eq.(p5|) is satisfied. 

• The inverse the matrix IS 3 {uJij, kij) must exist. 

From this result we can conclude that provided a solution (|3(]) exits, the integrability of the NPDE is completely 
determined by the linear part of the evolution equation. These are also the conditions which guarantee the integrability 
of Eq. (p8[) . In the following section, we apply these concepts to examine the integrability of some NPDEs. 



4 Examples 

In this section, we will apply the machinery developed in the precious sections to investigate the integrability of various 
NPDEs. As an first example, we consider the nonlinear Schrodinger equation: 

idtu = d xx u + 2uu*u (47) 

If we substitute: 

00 

u(x, t) = e i<™ e i(,* 2 -b 2 )t e i<t> A n e- n(bx - 2abt ^ (48) 

71=1 

into Eq.(ff7j), we obtain: 

oc oc n— 2 n— I— 1 

[(1 - n 2 )b 2 ] A n e- n ^- 2abt ^ = 2 ]T ]T £ A^A^^-^- 2 ^ (49) 

ri— 1 n — 1 1—2 rn — 1 
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It can be verified that for n = 1 the linear dispersion relationship u> = —k 2 (k = a + hi) is satisfied. Since both the 
left-hand side and the right-hand-side can be expanded in the same Fourier basis-functions, we can determine the 
expansion coefficients by the following recursion relationship: 

L (n) (k)A n = R {n) ; k = a + bi (50) 

where 

n — 2 n—l—1 

LW = [1 - n 2 ]b 2 - RM =2^2 J2 AiA m A n - m _i; (51) 

1=2 m=l 

If we assume that A\ — A, than by computing all the coefficients A n , and carrying out the summation, similarly as in 
Eq.(Ttl]), we obtain the NLS-soliton: 



i(x, t) = 2be lax e l{a2 - b2 ' >t e l ' t >sech(bx - 2abt + £ ) £o = log 



(52) 



Similarly as for the KdV equation the two-soliton solution of the nonlinear Schrodinger equation can be computed by 
considering solutions: 

oo 

u{x,t) = e ^+^ e i{al+al-bl-bl)t e i4> C{n,m)e~ n{[bl+b2]x - 2[aibl+a2b2]t) (53) 

By generalizing this procedure as presented in Sec. 3, the N-soliton solution of the nonlinear Schrodinger equation can 
be computed. 

As a second example we consider the coupled nonlinear Schrodinger equation: 



iuu = u lxx + {\ui\ 2 + \u 2 \ 2 )ui = 
iu 2t = u 2 xx + (\u2\ 2 + \ui\ 2 )u 2 = 



2 , |.,.|JU._„ (54) 



If we make the following substitution for the solution u(x,t) = [ui(x, t), u<i{x, t)] T : 

oo 

u(x, t) = e ^ e ^ 2 ~b 2 )t J- An ^n(^2aU) A („) = ( A W )A V)y (55) 



into Eq.(54), it can be easily verified that both the left-hand side and the right-hand side of Eq. (|54|) can be expanded 
in the same basis-functions. This is due to the fact that both u\(x, i) and u 2 {x, t) have the same dispersion relation 
ui(k) — ~k 2 . As a result, we can determine the expansion coefficients A^") by solving the following recursion relation 

L («) (fc) A<"> = R (n) k = a + bi (56) 

where 

B-an-J-l / Al)Am) An-m-l) Al) Am) An-m-l) \ 

t (") T T 1 n 2 l h 2 - R ™ — 1 1 1 + 2 2 1 • ^71 

L -l[L-n\b, K -2^2^ Al) Am) An-m-l) , Al) Am) An-m-l) > V° ' ) 

1=2 m=l V A l A l ^2 +^2 A 2 A 2 J 

As a more general example, we consider the Manakov problem which plays a role in optical fiber communication ^ |^] : 

iu u + iSui x = u lxx + (|ui| 2 + \u 2 \ 2 )ui = , , 

iu 2t - iSu 2x = u 2xx + {\u 2 \ 2 + K| 2 )u2 = 
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The Manakov equation (|5q ) reduces to the nonlinear Schrodinger equation if 5 = 0. For a non-negative S however, u\ 
and u 2 have different linear dispersion relations. If we substitute: 



u(x t) — e i ( a i+ a 2) x e ~ i ( a i~ a 2) St e i ( a i+ a 2~ b i~ b 2) t e il t> \ * C(n m s j e ~ nb i( x ^ St + 2a i t ) e ^ mb 2( x + St + 2a ^ t ) 



n,m— 1 



(59) 



into Eq.(47), we obtain that the linear dispersion relationships are determined if the coefficients Ci(l, 0) and 6*2(0, 1) 
are unequal to zero. It can be verified by substitution that in this case the following linear dispersion relations are 
valid: 

Wifa) = 8k, - k 2 k t = a, + ib, i £ {1, 2} (60) 

In general the coefficients C(n, m) are determined by the following linear algebraic relationship: 

~L ( - n > m \k 1 ,k 2 )C(n,m) = R,(»> m ) (61) 

where 

L (n ' m) (/c!, k 2 ) = (1 - n 2 )b\ + (1 - rr?)b\ + 2i {ma x b 2 + na 2 b l T m(5 2 8} + 2a x a 2 - 2nmb 1 b 2 T 28a 2 (62) 

and 



n — 2 n —2n—l—ln—l—l 

R (n ' n0 = EE E E 

l =2 l' =2 m=l m' = l 



Ci (I, l')Ci [m, m')C x (n-m-l,ri -m' - V) 
d (I, l')Ci (to, m')C 2 (n -m-l,ri - m' - V) 

C 2 (l, l')C 2 (m, m')Ci(n - to - l,ri - m! - V) 
C 2 ll,l')C 2 (m,m')C 2 ln - m - l,n' - m! - V) 



(63) 



In principle general solutions for the Manakov problem can be obtained by following the procedure as described in 
Section 3. It is interesting to conclude that although the Manakov problem has no soliton solutions, the existence of 
the solutions (^9|) seem to suggest that the the Manakov problem is integrable. The discussing about the integrability 
can be made general in a similar manner as presented in Section 2. The most simple solution of the Manakov equation, 
however has four spectral parameters. This makes that this equation has no soliton solutions, but the most simple 
solution is an interacting 2-soliton solution. From a physical point of view, the Manakov problem describes the 
propagation of optical solitons over their principal birefringence axis. It is well know that optical solitons propagating 
through a birefringent medium are unstable. This can be explained by the fact that initially localized optical pulses 
disintegrate due to their soliton soliton interaction. 

As a last example, we consider the three dimensional nonlinear Schrodinger equation: 



idtu 



3 

^2 d 2 Xn u + 2uu*u 



(64) 



If we substitute: 



U(x,t) = e ^-x e i(a.a-b-b)t e i0 j4ne -n(b-x-2a.bt) 



(65) 



into Eq. fl47|) , we obtain: 



00 

E 

71=1 



00 n — 2 n — l — 1 



[(1 - n 2 )b ■ b] A n e-^ h — 2 ^ = 2J2J2J2 AA^^e-^— 2 ^ 



(66) 



n—1 1—2 m—1 
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In Eq.(|35|) and Eq.(|66]), it is used that x = (xt, X2, X3) 71 , a = (at, a^, az) T and b = (pi, £>2, bs) T . It can be verified that 
for n = 1 the linear dispersion relationship lo 2 = — k • k(k = a + hi) is satisfied. Since both the left-hand side and the 
right-hand-side can be expanded in the same Fourier basis-functions, we can determine the expansion coefficients by 
the following recursion relationship: 

L< n > (k) A n = R (n) ; k = a + hi (67) 

where 

n— 2 n-l—l 

7>)(k) = [l-n 2 ]b-b; RW=2J2 E A l A m A n-m-l\ (68) 

1=1 m=l 

Similarly as in the one-dimensional case, explicit solutions of the three dimensional nonlinear Schrodinger equation can 
be obtained by carrying out the summation of the expansion coefficients. The discussion can be made more general 
by using more general expansion functions, similarly as in Eq.(|53|) 

5 Conclusions 

We have presented a novel method to test the integrability of large classes of NPDEs. We can conclude that the 
integrability of NPDEs depends on two important conditions. The first condition is that the NPDE must have a 
solution of the form (^TJ), or more precisely, the nonlinearity must have the property ([31]). This ensures that the 
nonlinear part of the NPDE can be expanded in the same basis functions as the linear part. If this is the case, the 
expansion coefficients of the solution can be determined by a simple linear algebraic equation. 

The machinery developed in this paper is derived by discussing the existence of N-soliton solutions of the KdV 
equation as an example. These results have lead to an insight which is used to present a general discussion with respect 
to the existence of N-soliton solutions of the KdV equation. In Section three the machinery is firstly generalized for 
NPDEs which have solutions that can be expanded in Fourier basis-functions and later generalized to general basis- 
functions. Secondly, the power of the method is demonstrated to derive solutions of the nonlinear Schrodinger equation 
in a simple way. It is shown that we can systematically derive large classes of solutions of the nonlinear Schrodinger 
equation. We have also demonstrated that it is possible to systematically derive solutions of the coupled nonlinear 
Schrodinger equation. Finally we derive solutions of the Manakov problem. The Manakov problem cannot be solved 
by inverse scattering techniques, however we have demonstrated that large classes of solutions can be derived. It 
is interesting to remark that the most simple solution of the Manakov problem involves two independent complex 
spectral parameters. This suggests that the most simple solution consists of two interacting solitons. The latter is 
related to the well known instabilities of the solutions of the Manakov problem. Finally it is indicated the the methods 
presented in this paper can be used to find solutions and to investigate the integrability to higher-dimensional NPDEs. 

As a concluding remark we want to state that the linearization principle as presented in this paper can be applied 
to large classes of NPDEs having a polynomial type of nonlinearity. It is well possible that this concept can be 
generalized to much larger classes of nonlinear problems. As an example we mention Rcf. Q, in which it is shown that 
a similar approach can be used to solve inverse scattering problems. 
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